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LONG-TERM  GOALS 

The  long-tenn  scientific  goal  of  my  research  is  to  better  understand  the  distribution  of  phytoplankton 
in  the  world’s  oceans  through  remote  sensing  their  influence  on  the  optical  properties  of  the  water. 

OBJECTIVES 

The  long-term  scientific  goal  of  my  research  has  been  focused  on  the  remote  sensing  of  phytoplankton 
and  other  constituents  from  measurements  water-leaving  radiance  propagated  to  the  top  of  the 
atmosphere.  Optically,  phytoplankton  and  other  constituents  reveal  their  presence  through  their 
influence  on  the  inherent  optical  properties  (IOP’s)  of  the  water.  The  main  effect  of  phytoplankton  is 
to  increase  the  absorption  of  light  by  virtue  of  the  strong  absorption  by  their  photosynthetic 
(chlorophyll  a)  and  accessory  pigments.  Other  constituents,  e.g.,  coccolithophores  reveal  theirs 
through  high  backscattering  (large  water-leaving  radiance).  Although  techniques  for  measuring  the 
absorption  coefficient  directly  (e.g.,  in-situ  AC9  measurements  or  in-vitro  filter  pad  absorption)  have 
become  accepted  by  the  scientific  community,  laboratory  techniques  for  measuring  backscattering  are 
tedious  and  subject  to  error,  and  in-situ  techniques  for  backscattering  are  in  their  infancy.  Thus,  there 
has  been  considerable  effort  devoted  toward  indirectly  inferring  these  IOP’s  by  virtue  of  their  influence 
on  the  apparent  optical  properties  (AOP’s),  e.g.,  the  diffuse  reflectance  of  the  water  (the  color  of  the 
water)  or  the  downwelling  irradiance  attenuation  coefficient.  These  AOP’s  are  perhaps  the  most 
frequently  measured  quantities  in  hydrologic  optics.  Clearly,  interpretation  of  such  observations 
requires  a  detailed  understanding  of  the  influence  of  phytoplankton  on  the  IOP’s  and  their  link  to  the 
AOP’s. 

Although  the  research  reported  here  is  centered  on  deriving  the  IOP’s  from  measurements  of  the 
AOP’s,  which  is  important  in  that  IOP’s  determined  from  AOP’s  are,  by  definition,  sampled  at  a  scale 
appropriate  for  radiative  transfer  and  for  remote  sensing.  A  secondary  focus  has  been  to  try  to 
understand  the  influence  of  particle  shape  on  light  scattering  by  marine  suspensions.  For  this,  we 
chose  to  attempt  modeling  the  scattering  of  detached  coccoliths  of  E.  huxleyi. 


APPROACH 

In  the  mid  1990’s  the  PI  and  coworkers  developed  techniques  for  retrieving  the  optical  properties 
(single  scattering  albedo  and  scattering  phase  function)  of  aerosols  from  measurements  of  sky  radiance 
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over  the  oceans  [Wang  and  Gordon,  1993;  Gordon  and  Zhang,  1995;  Zhang  and  Gordon,  1997a, 
1997b;  Cattrall,  Carder,  and  Gordon,  2003].  These  were  the  basis  for  development  of  algorithms  for 
inverting  irradiance  and/or  radiance  profiles  in  the  water  to  obtain  profiles  of  the  absorption 
coefficient,  a,  and  the  backscatterinng  coefficient,  b/,.  The  algorithms  are  iterative  in  nature  with  the 
radiative  transfer  equation  (RTE)  solved  during  each  iteration  using  trial  IOP’s.  The  key  was 
conceiving  a  simple  method  of  improving  the  estimated  IOP’s  after  each  iteration  by  comparing  the 
predicted  and  measured  AOP’s.  When  the  latter  agree  within  experimental  error,  the  IOP’s  were 
deemed  to  be  those  of  the  medium. 

For  the  study  of  light  scattering  by  irregularly  shaped  particles,  the  discrete-dipole  approximation 
[Draine  and  Flatau,  1994]  was  used,  and  the  coccoliths  were  modeled  as  disk-like  particles. 

WORK  COMPLETED 

Table  1  provides  a  summary  of  the  inverse  problems  that  have  been  examined  during  the  course  of  this 
study.  These  are  detailed  in  a  series  of  papers.  The  first  paper  reported  development  of  inversions  for 
homogeneous  media  with  elastic  scattering  [ Gordon  and  Boynton,  1997],  the  second  extended  the 
development  to  elastically-scattering  vertically-stratified  water  bodies  [ Gordon  and  Boynton,  1998], 
the  third  treated  inversion  in  the  presence  of  the  interfering  effects  of  Raman  scattering  into  the 
spectral  band  of  interest  [Boynton  and  Gordon,  2000],  and  the  fourth,  provided  an  improvement  to  the 
Gordon  and  Boynton  [1998]  algorithm  for  clear  water  (i.e.,  waters  in  which  scattering  by  the  water 
itself  makes  a  significant  contribution  to  bfz),  where  z  is  depth)  [Boynton  and  Gordon,  2002].  For  a 
complete  program  of  AOP  inversion  in  natural  waters,  the  only  process  we  have  overlooked  is 
fluorescence,  i.e.,  inelastic  scattering  into  a  given  band  of  wavelengths  from  all  shorter  wavelengths. 

A  method  for  solving  this  problem  was  provided  by  Gordon  [2001]  and  applied  to  simulated  data  by 
Boynton  and  Gordon  [2002]. 

Table  1.  Matrix  of  Inverse  Problems  Relevant  to  Hydrologic  Optics 


Problem 

Water  Body 

Homogeneous 

Stratified 

Elastic  (Depth  =  oo) 

Z 

Z 

Elastic  (Depth  <  oo) 

Z 

© 

Raman  (Depth  =  oo) 

Z 

Z 

Fluorescence  (Depth  =  oo) 

Z 

X 

Z ;  Inverse  solution  completed 

*:  Inverse  solution  not  likely  to  be  possible 

©:  Planned  approach  failed 
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A  theoretical  study  modeling  light  scattering  by  detached  coccoliths  as  disk-like  structures  was  carried 
out  by  Gordon  and  Du  [2001], 

Finally,  as  a  confluence  of  our  work  on  inverse  methods  and  coccolith  scattering,  during  the  current 
reporting  period  we  used  the  inversion  methods  that  we  developed  to  retrieve  b/,(z)  as  a  function  of 
wavelength  and  coccolith  concentration  by  inverting  irradiance  data  taken  in  intense  coccolithophore 
blooms.  A  paper  describing  this  effort  is  in  preparation. 

RESULTS 

The  entries  in  Table  1  summarize  the  results  of  our  study  on  inverse  methods  in  hydrologic  optics. 
Typically,  the  inversions  provide  IOP’s  that  reproduce  the  input  AOP’s  almost  exactly  or  ,  with  field 
data,  within  experimental  error.  The  check  marks  indicate  the  completion  of  an  inversion  method  and, 
except  for  the  case  of  fluorescence,  a  successful  application  to  field  data.  The  inversion  technique  we 
developed  for  fluorescence  [ Gordon ,  2002,  Boynton  and  Gordon,  2002],  although  functional  with 
synthetic  data,  is  very  unstable  as  inversion  in  the  presence  of  fluorescence  is  ill-conditioned.  The 
problem  can  be  solved  (at  least  by  our  techniques)  only  when  the  fluorescence  is  strong  —  a  situation 
not  encountered  in  the  natural  environment.  Initially,  we  believed  that  we  had  found  a  viable  approach 
to  inversions  in  stratified  waters  in  which  the  bottom  made  a  significant  contribution  to  the  upwelling 
light  field.  Unfortunately,  we  were  unable  to  separate  bottom  reflection  from  backscattering,  e.g.,  from 
a  layer  of  particles  just  above  the  bottom,  and  this  problem  remains  unsolved. 

Our  inversion  of  b/,(z)  as  a  function  of  wavelength  and  coccolith  concentration  in  an  intense 
coccolithophore  bloom  yielded  a  backscattering  spectrum  that  was  similar  to  one  of  the  Gordon  and 
Du  [2001]  models  —  most  complex  shape  studied  (the  "fishing  reel").  Figure  1  shows  an  SEM  image 
of  a  detached  coccolith  (left  panel)  along  with  the  “fishing  reel”  model  (right  panel),  which  can  only 
approximate  the  gross  structure  of  the  particle.  Figure  2  compares  the  spectral  variation  of  the  retrieved 
backscattering  coefficient  in  the  surface  layer  (red  lines)  with  that  predicted  for  the  “fishing  reel” 
model.  The  most  interesting  aspect  of  the  comparison  is  the  increase  in  scattering  at  wavelengths 
greater  than  550-600  nm  observed  in  both  the  model  and  the  measurements.  This  increase  in 
backscattering  is  due  to  constructive  interference  between  rays  internally  reflected  by  the  two  plates 
(resonant  scattering).  Finally,  although  the  spectral  variation  of  bb  is  well  reproduced,  the  absolute 
magnitude  of  the  backscattering  on  a  per  coccolith  basis  is  smaller  (~  factor  of  2)  than  observations, 
and  on  a  per  calcite  basis  is  somewhat  higher  than  field  measurements. 

IMPACT/APPLICATIONS 

The  AOP  inversion  algorithms  that  were  developed  under  this  grant  provide  realistic  profiles  of  a  and 
The  research  also  shows  that  although  backscattering  by  disk-shaped  objects  is  very  dependent  on 
morphology,  it  can  be  understood,  and  that  models  of  irregularly  shaped  particles  for  light  scattering 
by  marine  particles  show  some  promise  for  understanding  scattering  properties  of  suspended  particles. 
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Figure  1.  SEM  image  of  a  detached  coccolith  from  E.  Huxleyi  (left  panel).  Fishing  reel  model  of 
the  coccolith  resembling  two  washers  of  outside  diameter  (2.7  (am  )  separated  by  a  small  gap 

(0.3  |am). 


A  (nm) 


Figure  2.  Comparison  of  the  spectral  variation  of  the  retrieved  backscattering  coefficient  in  the 
surface  layer  (red  lines)  with  that  predicted  for  the  Gordon  and  Du  [2001]  “fishing  reel ”  model. 
The  most  significant  aspect  of  the  comparison  is  the  increase  in  scattering  at  wavelengths  greater 
than  550-600  nm  observed  in  both  the  model  and  the  measurements 
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RELATED  PROJECTS 


We  are  collaborating  with  A.  Morel  to  apply  our  irradiance  inversion  algorithm  to  analysis  of  the 
OLIPAC  data  set  and  with  the  T  J.  Smyth  (PML)  in  applying  our  algorithms  to  data  acquired  in 
coccolithophore  blooms  to  further  delineate  their  optical  properties.  We  have  also  been  collaborating 
with  R.  Leathers  (NRL)  in  comparing  retrieval  algorithms  operating  on  irradiance  data  from  the  Gulf 
of  California. 
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